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Abstract 

We investigate the coherent manipulation of interacting Rydberg atoms placed inside a high- 
finesse optical cavity for the deterministic preparation of strongly coupled light-matter systems. We 
consider a four-level diamond scheme with one common Rydberg level for N interacting atoms. One 
side of the diamond is used to excite the atoms into a collective ‘superatom’ Rydberg state using 
either 7r-pulses or stimulated Raman adiabatic passage (STIRAP) pulses. The upper transition on 
the other side of the diamond is used to transfer the collective state to one that is coupled to a field 
mode of an optical cavity. Due to the strong interaction between the atoms in the Rydberg level, 
the Rydberg blockade mechanism plays a key role in the deterministic quantum state synthesis of 
the atoms in the cavity. We use numerical simulation to show that non-classical states of light 
can be generated and that the state that is coupled to the cavity field is a collective one. We 
also investigate how different decay mechanisms affect this interacting many-body system. We 
also analyze our system in the case of two Rydberg excitations within the blockade volume. The 
simulations are carried out with parameters corresponding to realizable high-finesse optical cavities 
and alkali atoms like rubidium. 

PACS numbers: 32.80.Ee, 42.50.Ex, 37.30.+i, 42.50. 
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I. INTRODUCTION 


We investigate using Rydberg atom interactions to deterministically synthesize collective 
quantum states that can be strongly coupled to an optical cavity for applications in quantum 
information science, quantum optics and the production of quantum light fields. The method 
enables a single or several collective excitations to be prepared in a high-finesse optical 
cavity The number of excitations is precisely controlled so that photon number states can 
be produced using the coupled atoms-cavity system. The system we describe in this paper 
differs from typical Jaynes-Cummings (JC) and Tavis-Cummings (TC) physics in that strong 
atom-atom interactions play a key role in the initial state preparation, which is the focus of 
this paper. 

Rydberg atoms have many exotic properties. For example, Rydberg atoms are large, 
interact strongly with one another, have relatively long lifetimes, and can be manipulated 
with external electric and magnetic fields [1]. Particularly important to the work presented 
in this paper, the interactions between ultracold Rydberg atoms can be larger than typical 
frequency stabilized laser spectral band widths and kinetic energies at distances of over 10 /um 
[2, 3]. The interactions can be manipulated easily with electric fields, relative to valence 
states [4]. There is a broad range of behavior that can be induced when the multi-level 
character of the interactions is taken into account. 

Rydberg atom blockade is a result of the strong interactions between Rydberg atoms. 
In Rydberg atom blockade, the long range interactions between Rydberg atoms suppress 
multiple excitations in a volume in which the interactions are large enough to shift the 
excitation laser out of resonance with the Rydberg atom transition frequency [5, 6]. It 
has been realized over the last 15 years that Rydberg atom blockade is useful for a host 
of applications, particularly in quantum optics and quantum information science. Rydberg 
atom blockade can be used to control quantum dynamics and prepare interesting quantum 
states [7-9]. Rabi oscillations have been observed in effective two-level ultracold Rydberg 
atoms [10-15] demonstrating long coherence times and collective behavior. Collective Rabi 
oscillations on Rydberg transitions have even been observed in a thermal rubidium (Rb) 
vapor on a time scale below 1 ns [16]. 

Atoms contained in cavities, both optical and microwave, have been important systems 
for investigating the quantum-classical boundary in physics. The ability to entangle atoms 
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with light presented by cavity quantum electrodynamics has also served as a model system 
for applications in quantum information science [17, 18]. A cavity quantum electrodynam¬ 
ical system with a single excitation or controlled number of excitations is useful both for 
fundamental studies [19] and applications that use quantum behavior, particularly those 
involving quantum entanglement [20-22], Entangled states can be created between distant 
atoms in one cavity or separate cavities [23-27]. Trapping of a single atom in a high-finesse 
optical cavity has been achieved [28, 29]. However, it is technically challenging to couple 
the light in the cavity to a single atom because of the small absorption cross-section [30-33]. 
The problem of small absorption cross-sections can be handled by using an atomic ensemble, 
where light more strongly interacts with collective atomic states. 

In this paper, we theoretically show that one can excite collective N atom superposition 
states using Rydberg atom interactions that can be strongly coupled to a high-finesse cavity. 
We observe y/N enhancement of the Rabi oscillations within the cavity showing that the 
collective state can be realized [34, 35] and analyze the decay mechanisms that affect this 
interacting many-body system. Collectively enhanced cooperativity factors help to improve 
the efficiency of photon generation out of the cavity. The efficiency can be improved by 
increasing the number of atoms. We also show that by chirping the laser excitation pulses 
in time, the single atom Rydberg blockade effect can be suppressed and we can create two 
Rydberg excitations within the blockade radius that are coupled to the cavity. 

We use a four-level diamond type atomic scheme as shown in Fig. 1(b). One side of 
the diamond is used to collectively excite the atoms to the Rydberg level using a pair of 
TT-pulses or a pair of counter-intuitive stimulated Raman adiabatic passage (STIRAP) pulses 
with Rabi frequencies Pi(t) and P 2 (t). By using a laser pulse, and the single atom 

cavity coupling strength, g, the state of the collective atoms-cavity system can be controlled. 
Rydberg atom blockade plays an important role in this process, since we are interested 
in cases where a single or controlled finite number of coherently shared excitation (s) are 
produced in a collection of 5 P 3/2 state Rb atoms. The N interacting atoms-cavity system 
is within the reach of experimental realization [36, 37]. The system can be used to prepare 
single photons using a single excitation and more complex quantum fight fields by preparing 
multiple excitations in the cavity. 

Recently, a JC model in the optical domain has been proposed for investigating Rydberg 
blockaded atomic ensembles in a cavity. The transmission through the cavity was studied 
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numerically by Monte-Carlo simulation [38]. In contrast to this work and others, our studies 
resemble a pumping scheme for inverting the population of a laser, except the upper state 
of the lasing transition is prepared in a collective state. The pumping scheme is similar to a 
Raman or superradiant laser and work done on low atom number masers and micromasers. 
However, for these types of lasers and masers, the atoms are coherent but non-interacting 
[39-41]. Our approach is useful for single photon generation [42, 43], as well as the genera¬ 
tion of other types of quantum light fields, and coherent optical manipulation for quantum 
information processing [25, 26, 44-47]. The collective atoms-cavity system can also be im¬ 
portant for generating unique quantum states such as superposition states of the atomic 
ensemble [48, 49]. The approach does not require phase matching and it may be easier to 
use this system than those based on phase matching to create more elaborate quantum light 
fields [50, 51]. Creating arbitrary quantum light fields remains an important challenge for 
quantum optics, so it is useful to explore a wide range of different strategies. The method 
that we explore in this paper may have technological and experimental advantages because 
the collective state is prepared on a set of transitions which have significantly different ex¬ 
citation wavelengths, > 15 nm for all transitions in Rb, than the cavity mode to which it is 
coupled. The cavity output can be filtered and detected with low background. 


II. MODEL 

Fig. 1 shows a conceptual drawing of the experimental setup. The setup consists of a 
collection of atoms trapped in a high-finesse optical cavity. The atomic energy level scheme 
is a four-level diamond scheme with atomic states denoted by \g), |e), |s) and |r). hwi is the 
energy of the atomic level, therefore the transition frequencies Uij — Uj—Ui where i and j 
are labels corresponding to the respective states \g), |e), |s) and |r). The Rabi frequencies are 
taken to be time dependent since we are concerned with applying pulses for the preparation 
of a collective state that we can couple to the cavity as well as with controlling the photons 
emitted from the cavity. The decay processes that we consider explicitly in our model are 
also shown. We envision the atoms to be held in the cavity with a far off resonance dipole 
trap (FORT) or a magnetic trap. The dipole trapping wavelength can ideally be set to a 
magic wavelength so that the frequency of the \g) -H- |e) transition is fixed in the cavity. 
We discuss example design parameters in this section as well as the following 2. Additional 
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decay processes and the consequences of the main approximations used are also discussed 
later in the paper. 


(a) Rydberg atoms-cavity system (b) Atomie-level scheme 



FIG. 1. (a) A Rydberg atom collective state in an optical cavity indicating Rydberg blockade, 
(b) Level scheme of a four-level atomic system driven by three laser pulses on transitions |gf) •(-)■ 
|s), |s) •(-)■ |r) and |r) •(-)■ |e) with Rabi frequencies Pi{t), P 2 {t) and respectively. The transition 
\g) -H- |e) is conpled to an optical cavity with coupling strength g. F^ is the decay rate of the 
Rydberg state, F^ = Fo/2 is the transverse spontaneons decay rate of the |e) state and k is the 
cavity decay rate. Energy levels are shown for Rb. 

One side of the diamond is used to collectively excite the atoms to a Rydberg level 
using a pair of vr-pulses or counter-intuitive STIRAP pulses with Rabi frequencies Piit) and 
P 2 {t)- Pi{t) drives the \g) -H- |s) transition, frequency cugs, while P 2 {t) drives the |s) -H- |r) 
transition, frequency cUgr- If the interaction between the Rydberg atoms is strong enough so 
that the \g) -H- |r) transition is shifted out of resonance after the first Rydberg excitation 
occurs, further population of |r) is blocked. 

The other side of the diamond system is used to produce a collective state that is close 
to resonance with a field mode of a high-finesse optical cavity. The optical cavity mode 
is coupled to the \g) -f-)- |e) transition since the resonant frequency of the cavity, cUc, is 
assumed to be close to the atomic transition frequency, ojge- In fact, cUc = ojge for most of the 
calculations that follow. The control laser is treated classically and has a frequency ui that 
couples to the |e) -H- |r) transition at frequency ui^r with a Rabi frequency ri{t). The cavity 
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field is treated as a quantum field. The single atom states of the coupled atoms-cavity field 
system are products of |e) or \g), each coupled to a ladder of photon number states. 

The density matrix equations are not easily solvable for an arbitrary number of atoms 
because the number of states increases dramatically with N. One way to reduce the basis set 
size in a many-body system is to use a mean-field approximation. For an ultracold Rydberg 
gas this means each highly excited atom shifts the Rydberg levels of surrounding atoms out 
of resonance with the laser field [52]. In this approach, atomic coherence is damped out 
and the quantum master equations reduce to many-body rate equations. This method is 
useful for understanding excitation suppression due to Rydberg atom blockade in cold gases 
[53]. A second way to practically address the increasing basis set size is to utilize the idea 
of local blockade and truncate the basis set when the number of excitations in the system 
exceeds some limit [54, 55]. The second method is appropriate when the system dynamics 
are constrained to a subspace of the complete basis. Rydberg atom dipole blockade naturally 
justifies this approach. In the spirit of the latter approach, we carry-out our calculations in 
a truncated Hilbert space, similar to that used by Saffman and Walker [56]. 

To prepare the initial collective excitation of the atoms, Pi{t) and P 2 {t) are detuned from 
the atomic transition as shown in Fig. 1(b). Ag is the detuning from the intermediate level 
|s), assumed to be fixed in time, and A(t) is the detuning from the Rydberg level |r). The 
one photon detuning A^ is taken to be larger than either of the excitation Rabi frequencies 
for excitation by vr-pulses, Ag ^ Ug — oJg — > Pi(t), P 2 {t) and the detuning satisfies the 

two photon condition, 

^P2 T ^Pl (^) 

where are the laser frequencies associated with Pi{t) and P 2 {t), respectively. Under 

these circumstances, the residence time in level |s) is small and |s) can be adiabatically 
eliminated. For STIRAP pulses, the one photon detuning A* has a very small effect and 
the required conditions are a 2-photon resonance between initial and final states and adia¬ 
batic evolution of the pulses. The adiabatic approximation means that the Rydberg state 
excitation can be replaced by a single photon excitation in this setup. However, most ul¬ 
tracold Rydberg work uses a 2-photon scheme like this one and a 2-photon excitation from 
the ground state of an alkali atom allows access to nS states, which generally have simpler 
interactions, or nD states, which have more sub-states to use. 

The effective Hamiltonian for the excitation side of the diamond system can be written 
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in a collective state basis as 


i/i = 5 [-2(A(t) + A.)|fi)(fi| - 4(A(t) + A.)|fiR>{RR| + {Si(()|G>{fi| + 

+H.C.}], (2) 


where 


Slit) 


y/NSit) 


S2it) 


VN^Sjt) 
1/2(2 +Ait)/A,) 


and Sit) 


Plit)P2it)/As, 


IG*) = b) = \9i-9i-9N), 

N N 

\R) = 1/vW \r/ = l/y/NY. W-P-gN), 

1=1 1=1 


( 3 ) 

( 4 ) 

( 5 ) 


N N 

\RR) = l/i/Pp = V\/^ \9i---ri---rj---gN), (6) 

+j=i +j=i 

and 

N\ 

^ (N - 2) \ 2 V 

The Hamiltonian for the emission side of the diamond system for N interacting atoms with 
truncated Hilbert space bounded by the two Rydberg excitations is 


H2 = h/2 2A,(|E,0)(E,0| + |E, 1)(E, 1| + 2|EE, 0)(EE, 0| + |ER, 0)(ER, 0|) 


+ 


R, 0)(E,0| + |R,1)(E,1| + y2|EE,0)(ER,0| + y2|ER,0)(RR,0 
+g (\/]V|G, 1)(E,0| + \/]V|G,2)(E, 1| + V(iV-l)|R, 1)(ER,0| 


V2(A^-1)|E,1)(EE,0 


H.c. 


( 8 ) 


where flit) is the Rabi frequency of the control laser, C^C) is the creation (annihilation) 
operator for the near resonant cavity mode, Ac is the cavity detuning, g is the coupling 
strength of one atom to the cavity and 

N N 

|E) = l/^Y. Ig) = (9) 

1=1 1=1 


N 

|EE) = l/y/]^ ^ |e,e,), 
+j=i 


( 10 ) 
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and 


N 

|ER) = l/^/m (11) 

i^j=l 

On the emission side of the system, the first element of the ket represents the collective 
atomic state. The second number of the ket represents the Fock, or photon-number state, 
of the near resonance mode of the cavity. |E) and |EE) are the collective states we are 
interested in preparing in the cavity [57]. Although much of the interesting physics in this 
system is due to the cavity quantum electrodynamical system, it is important to note that 
we are also addressing the preparation of the atomic collective state and the effect of several 
different decay mechanisms. As a result, both Hi and H 2 are important. Stark shifts caused 
by the excitation lasers can be compensated with detunings. 

In our calculations, we assume that the interaction volume is less than the blockade 
volume. This allows for the modeling of the blockade effect through a straightforward 
detuning of states with more than a single Rydberg excitation. The level shift between a 
pair of atoms labeled i and j each in a Rydberg state can be estimated using 

A 

where two types of simplified interatomic interactions can be considered: (a) dipole-dipole 
interaction where p = 3, and (b) van der Waals interaction where p = 6. Cp are the dispersion 
coefiicients, and li and Ij are the positions of the and atoms, respectively. Here, 
we are primarily interested in repulsive potentials due to the atom-atom interaction. More 
realistic Rydberg atom interaction potentials have a complicated dependence on internuclear 
separation, including many avoided crossings that can be manipulated with electric and 
magnetic fields [2-4, 58, 59]. We use Eq. 12 to estimate the magnitude of the interaction 
whose exact value and radial dependence is not critical provided it is large enough for 
Rydberg blockade for the present discussion. The Hamiltonian for atom-atom interaction 
can be written as 

Ha,a = h/2'^\ri){ri\AR\rj){rj\. (13) 

i<j 

The condition that all atoms in the interaction volume he within the blockade volume 
translates to 

jnwt+lfaWP 

y2|Pi(()p + U/4 
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Ar » max 


(14) 



where is the spontaneous decay rate of the intermediate state |s). The situation where all 
atoms are within a blockade volume can be realized by using sufficiently narrow bandwidth 
lasers to form an interaction volume with Rydberg atom interactions strong enough to 
satisfy Eq. 14. In this work, the atoms are also assumed to be cold enough that they can be 
approximated as stationary, or frozen in space. These ideas are extendable to countable 
numbers of excitations by increasing the interaction volume or decreasing the blockade 
volume. We also consider the case where we want to compensate the shift due to the 
one Rydberg excitation blockade, Ar. This shift can be compensated by changing the time 
dependent detuning A(t) of the excitation fields, so that we can excite two Rydberg atoms in 
the blockade volume. More elaborate states can be created by, for example, applying more 
pulses to create several collective excitation states. Note that for 2 excitations, rather fiat 
avoided crossings that occur in Rydberg atom interactions can be used to create 2 excitation 
states with minimized spatial correlation [4, 58, 59]. 

The total Hamiltonian of the N atom system coupled to the cavity is 


HT = H^ + H2 + Ha,a. (15) 

The density matrix operator p of the N interacting atoms and all the fields obeys Liouville’s 
equation 

-^ =-^[Ht, p] +Ap, (16) 

where the decays are due to spontaneous emission from |e) and |r) and loss from the cavity. 
These effects are included in the second term of the equation using the Lindblad operator 
C. Cp can be written as 


Cp — Crp + T_lp + (17) 

with 

CrP = Vr{2UpLl-[ULlp]), (18) 

C^p = T^[2L^pLi - [L±4,p]), (19) 

, pj), (20) 
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Ll = \R,0){E,0\ + \R,1){E,1\ + V2\RR,0){ER,0\ + V2\ER,0){EE,0\ 


+ ^/2{N-l)\RR,0){R,0\, (21) 


L{ = ^N{N-1)\ER, 0)(L, 0| + ^N{N -1)/2\EE, 0)(L, 0| + y/N\E, 1){L, 0| 

+^/N\E,0){L,0l (22) 


and 


Li = |G',2)(G,1| + |G',1)(G',0| + \E,1){E,0\ + \R,1){R,0\, (23) 

where is the Rydberg state decay rate, F^ is the transverse spontaneous emission rate into 
the free space modes and k is the cavity decay rate. A dummy state \L, 0) is used to model 
decay into non-cavity modes and loss from the collective state manifold due to dephasing. 
Here the ‘0’ includes all other modes of the light field to simplify the notation. In particular, 
we use the dummy state to model the effect of additional dephasing mechanisms described 
in the Discussion section and to estimate their effect. 

We define the parameter Fq = 27 r x 6 MHz, the free space spontaneous decay rate corre¬ 
sponding to ®^Rb, as a scale for the decay rates used for the calculations that follow. The 
total Rydberg state decay rate is F^ ~ 27r x 1.4 kHz for n = 90. The black body radiation 
decay rate of the Rydberg state is very small, F;,;, ~ 27r x 0.2 kHz for n = 90, compared to 
the total decay rate. Consequently, we ignore black body decay which would cause atoms to 
leave our model system [60]. At large principal quantum numbers, n, where the interactions 
are strongest, which are most appropriate for the scheme described in this paper, the neglect 
of black body decay is best justified. F^ is small but we include it for completeness in the 
Liouville equations. 

The time variation of the electric field of the laser pulses is of the form, 

£;(t) = Ao sech[^^]. (24) 

T 

Aq is the amplitude, to is the time when the pulse reaches its maximum value and r is the 
pulse width. 27r x is taken to be smaller than the shift due to the interaction of Rydberg 
atoms Ar. 
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III. CAVITY AND INTERACTION PARAMETER ESTIMATES 


The key parameter for determining the emission rate into the cavity for a single atom is 
which is the ratio of the rate of emission into the cavity mode to the loss due to 
spontaneous emission into other electromagnetic held modes. This ratio is an expression for 
the Purcell factor, Fp = / (kT^), described in more detail in the next section, and is related 

to the single atom cooperativity factor, C = Fp/2. The cavity decay rate k is determined 
by the mirror rehectivities. The single atom coupling strength g is determined by the 
atomic transition dipole moment, /r, for the transition near the cavity resonance and the 
mode volume of the cavity, V[Lc, Rc,^]- Realistic cavity parameters limit Fp. For collective 
excitations of N atoms, the effective Purcell factor and collectively enhanced cooperativity 
increase, since p increases like \/]V. 

For our numerical calculations, we use parameters that have been achieved in recent 
experimental setups [8, 9, 15, 37] and choose ®^Rb to investigate the general system. To 
reach the strong coupling regime, we can achieve a single atom cavity coupling strength 
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UJr 


2heQV[Lc, Rci A] 

of up to 27r X 50 MHz for realistic cavity parameters. The cavity mode volume is 


(25) 


TT 




(26) 


and the beam waist is 


^o[Rc7 Rc'! 



(27) 


Rc is the radius of curvature of the mirrors, Lc is the cavity length and A is the wavelength of 
the cavity mode. We concentrate on numerical results corresponding to conditions near the 
strong-coupling regime, g » k, Fx,Fr as well as the leaky cavity regime, k » g^jK >> 
Fx,rr. Choosing realistic parameters, g and k, does not allow a clean separation of the 
strong-coupling and leaky cavity regimes. Relevant cavity parameters are length Lc ~ 50 
to 300 pm and beam waist wq 15 to 22 pm. A mirror radius of curvature of around 
Rc = 25 mm is suitable for the cavity. 

The blockade radius depends on the states of the Rydberg atoms. The blockade radius 
for principal quantum number n = 90, corresponding to the ^ 2% x 30 MHz excitation 
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Rabi frequency, is ~ 9 jim. At n = 90, the shifts due to the interaction of Rydberg atoms 
Aij ~ 27r X 220 MHz at 6.3 jim [2, 3]. For densities of ~ 10^^ cm“^, there are ~ 3 x 10^ atoms 
in the blockade volume with \/N shot to shot atom number fluctuations of ~ 50 atoms. The 
possible number of atoms in the blockade volume far exceeds the atom numbers considered 
in this paper. Small atom numbers, as studied here, can be generated by decreasing the 
trap density. 


IV. THE PURCELL REGIME 


The rate of spontaneous emission into the cavity mode can be enhanced by the interaction 
of the atoms with a matched resonant cavity, commonly called the Purcell effect. The Purcell 
regime occurs when k > g > V±. The enhancement, Fp, written in terms of the cavity mode 
volume and mirror reflectivities, Ri and R 2 , is [61], 

” V[L,,R„\]' ^ ’ 

where the cavity finesse is 


F{Ri,R2) 


27r 

(1 — Ri) + (1 — R 2 ) 


The decay rate into the cavity mode due to the Purcell effect is then 


(29) 


r, = Fpr„. (30) 

As a realizable example, we take Ri = 0.999985 and R 2 = 0.99985 with radius of curvature 
Rc = 25 mm. We can obtain a finesse F ~ 3.8 x 10^ at 780 nm with these reflectivities. For 
Purcell factor Fp ~ 10, the cavity length is Lc ~ 50 fim, which still allows the positioning of 
a 30 pm FORT inside the cavity. For this example, p 27r x 50 MHz and k ~ 27r x 72 MHz. 

In the Purcell regime, the coupling to the cavity, accounted for by g, can be replaced 
with an enhanced spontaneous emission rate into the cavity mode. We explore this regime 
by replacing the direct coupling of the atoms to the cavity with an incoherent decay term 
into one mode of the electromagnetic field 

Ll = ^/N\E,0){G,1\ + Vn\E, 1)(G,2| + ^/N^\ER,0){R, 1| 

+ V2{N-1)\EE,0){E,1\ (31) 
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with 


^pP — rp(2LppLj, — [I/pI/p, p]). (32) 

Effectively, the cavity is explicitly removed from the problem and the picture becomes one 
where an atom, or collective excitation, in free space preferentially emits into a specific mode 
of the electromagnetic field. The Purcell picture provides physical insight into preparing and 
detecting a countable number of excitations created in a cavity. 

V. VERIFICATION OF THE MODEL 

First, for low N and a single Rydberg excitation as a target, we consider the possibility 
of 2 Rydberg excitations occurring in our system. In these calculations, we show that, for 
Ar >> 0 and A = 0, there is insignificant leakage out of the Hilbert space bounded by a 
single Rydberg excitation. For calculations with N > 3, we use these results to justify the 
truncation of the Hilbert space at the 2 Rydberg excitation level. These results are also 
supported by numerical results of many other works done on Rydberg blockade. In these 
studies it is demonstrated that a large probability (> 0.98) for single Rydberg excitation 
of a superatom can be achieved while the probability of multiple excitations is negligible 
[46, 62, 63]. 

In Fig. 2, we show the results for N = 1 — 3 atoms with one Rydberg excitation. The 
model is efficient for investigating the case of V > 3, but here we focus our attention on 
a few interacting atoms within the cavity. The time dependence and peak amplitudes of 
the Rabi frequencies are shown in Fig. 2(a-c), while the probabilities for various states are 
shown in Fig. 2(d-i). For N = 1 and N = 2, the results are exact within the bounds of the 
model. To evaluate the solution for N > 3, we ignore the higher order excitation states. 
The population in those states are < 0.0015 justifying the truncation at the 2 excitation 
level for the work addressed in this paper. 

To further verify the model, we numerically investigate the preparation of a single col¬ 
lective excitation and the observation of a y/N enhancement of the Rabi oscillations of the 
superatom in the cavity. Collective states coupled to the cavity will exhibit a \/N enhance¬ 
ment in their effective atoms-cavity coupling strength, '/Ng. 

For V = 1, we used two tt pulses, Pi{t) and P 2 {t)i applied at the same time to prepare 
an atom in the Rydberg state. Subsequently, a control laser pulse with Rabi frequency Vt{t) 
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FIG. 2. (a-c) The time variations of the Rabi frequencies of the laser pulses applied to a four- 
level system for N interacting atoms. The effective Rabi frequency S{t) is represented by the 
black solid curve. The red dotted curve shows the Rabi frequency of the control laser pulse 
and the blue dashed curve shows the cavity coupling strength g. Panels (d-i) show the 
corresponding probabilities for population of states |G,0) (black solid), |i?, 0) (red dashed), \E^0) 
(blue solid), |G, 1) (magenta dot-dashed) and |T,0) (green dotted). Parameters for the calculation 
are P^ = 27r x 1.4 kHz, Pq = 27r x 6 MHz, k = 27t x 1.4 MHz = 0.1^, Ac = 0, Aji — 2 tix 220 MHz, 
A(t) = —27r X 110 MHz and As — 2 ti x 110 MHz. 
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is turned on to transfer the collective state to one that is near to resonance with the cavity. 
The population oscillates within the cavity between \E, 0) \G^ 1). Due to the cavity decay 

rate, k, = these oscillations decay with time [64], For N = 2 and N = 3, two counter¬ 
intuitive STIRAP pulses [62, 63] with pulse areas of tt were used to create a collective state 
with a single Rydberg excitation. We also increased the amplitude of the control laser pulse 
to optimize the transfer probability. In Fig. 2(e) and Fig. 2(f), the frequency of the Rabi 
oscillations are enhanced by '/N in the cavity. The decay rate and detuning of the cavity 
were the same as in Fig. 2(d). We also obtained the same results as shown in Fig. 2(e) and 
Fig. 2(f) using two superposed tt pulses (i.e. without using a STIRAP sequence). The plots 
in Fig. 2 demonstrate that only a single excitation is prepared in the cavity and it collectively 
interacts with the cavity held. This result shows our numerical approach is consistent with 
expectations. 


VI. NUMERICAL RESULTS 

In Fig. 2(g-i), we can see that the probability of populating the cavity ground state, 
IG, 0), is increased with increasing number of atoms. For N > 1, due to imperfect 7r-pulses, 
the population returning to the Rydberg state increases compared to Fig. 2(d-f) where 
there is no decay [45]. These results suggest that the probability of populating the cavity 
ground state depends on the number of atoms and the spontaneous emission rate out of 
the cavity. As we increase the coupling to the cavity mode, it, not surprisingly, makes the 
photon production from the cavity more efficient. This is the result of increasing the effective 
(collectively enhanced) cooperativity. 

Fig. 3 shows the time variation of the Rabi oscillations for N atoms in a resonant cavity, 
with Ac = 0 [Fig. 3(a-c)], and Ac = 0.3g [Fig. 3(d-f)], for different cavity decay rates 
K = 0.1,0.2,0.35'. As the cavity decay rate is increased, the Rabi oscillations are more 
quickly damped out [65]. However, increasing the number of atoms enhances the frequency 
of the oscillations by the a/]V factor as shown in Fig. 3(a-c). These oscillation frequencies 
are further enhanced by detuning the cavity [66] at the expense of a decreasing probability 
amplitude for photon production as shown in Fig. 3(d-f). These results are consistent with 
the JC model. 

The rate of single photon emission for a given number of atoms versus time is shown in 
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FIG. 3. Time variation of Rabi oscillations for N interacting atoms in the cavity for cavity detuning 
(a-c) Ac = 0 and (d-f) Ac = 27rx 4.3 MHz = 0.3^. Three different cavity decay rates are shown. 
The black solid curve corresponds to = 27rx 1.4 MHz = 0.1^, the red dotted curve corresponds 
to = 27rx 2.85 MHz = 0.2^, the blue dashed curve corresponds to = 27rx 4.3 MHz = 0.3^, and 
= 0. Other parameters are the same as in Fig. 2. 
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Fig. 4 with cavity decay rate k = O-l^" and Fj^ = ro/2. In this case, we observe oscillations 
in the rate of photon emission that change with N. The most interesting thing here is that 
the rate of photon emission is more strongly peaked at early times by increasing N. These 
results demonstrate that preparing the collective excitation helps to have the photon emitted 
at a specific time, due to an increase in the effective cooperativity from the collective state. 



FIG. 4. Rate of single photon emission efficiency for = 1 — 3 atoms with one Rydberg excitation 
vs time. The black solid curve represents the photon emission with = 1, the red dotted curve 
corresponds to A^ = 2 and the blue dashed curve corresponds to N — 3. The cavity decay rate is 
K = 2nx 1.4 MHz = t).lg and Fj^ = 27rx 3 MHz = Fo/2. Other parameters are the same as in 
Fig. 2. 

In Fig. 5, we show the results for N = 1—3 atoms for the full model and with an incoherent 
decay based on an enhanced spontaneous emission into the cavity mode as discussed in 
section IV. The time dependence and amplitudes of the Rabi frequencies for the incoherent 
decay associated with the Purcell effect are shown in Fig. 5(a-c), while the probabilities for 
various states are shown in Fig. 5(d-f). Rabi frequencies for the full model are shown in 
Fig. 5(g-i) with the associated probabilities for each state in Fig. 5(j-l). In both cases, we 
observe that the photon emission efficiency is enhanced with number of atoms because the 
cooperativity is effectively increasing due to the collectively enhanced /r. The efficiency, for 
a given number of atoms, for an incoherent decay for the Purcell effect based calculation is 
greater than in the full model because there is significant transient population of the cavity 
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FIG. 5. Panels (a-f) show the results with an incoherent decay in the Purcell regime. Panels (g-1) 
show the results with an enhanced spontaneous emission rate into the cavity mode in the weak 
cavity coupling regime, (a-c, g-i) The time variations of the Rabi frequencies of the laser pulses 
applied to a four-level system with N interacting atoms. The effective Rabi frequency S{t) is 
represented by the black solid curve. The red dotted curve shows the Rabi frequency of the control 
laser pulse and the blue dashed curve shows the cavity coupling strength g. Panels (d-f, j-1) 
show the corresponding probabilities for the population of states |G, 0) (black solid), \E^0) (red 
dot), |i?, 0) (blue dashed), |G, 1) (magenta dot-dashed) and |T,0) (green dotted). The cavity decay 

rate is = 27rx 72.6 MHz. Other parameters are the same as in Fig. 2. 
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excited state in the full calculation. Population of \E^ 0) can spontaneously decay into a free 
space electromagnetic field mode, degrading the probability that the photon is emitted into 
the cavity mode. 
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FIG. 6. Photon emission efficiency as a function of the cooperativity, C, for N interacting atoms 
with different cavity length, Lc- The mirror radius of curvature Rc = 10 cm. k is decreasing to 
increase C for fixed g. Other parameters are the same as in Fig. 2. 

In Fig. 6, the photon emission efficiency is plotted as a function of the cooperativity, C, for 
N interacting atoms with different cavity length Lc- As we decrease Lc, the photon emission 
efficiency can be improved because Lc increases the single atom g. The efficiency of photon 
emission also depends on the control laser pulse which transfers the Rydberg population 
into the cavity. For large cavity lengths, we can see that the photon emission efficiency 
first increases with C and then decreases. Although C is increasing, k is decreasing. The 
efficiency decreases because k becomes smaller compared to F^. As we increase the number 
of atoms, the photon emission efficiency can be further improved by increasing the effective 
cooperativity but only when cavity loss beats out spontaneous emission. 

In Fig. 7, we show the results for A = 1 — 3 atoms in the weak {k, > g > F^) and 
strong {g > K > T±) coupling regime. The time dependence and peak amplitudes of the 
Rabi frequencies for the weak coupling regime are shown in Fig. 7(a-c) and for the strong 
coupling regime in Fig. 7(g-i). The probabilities for various states are shown in Fig. 7(d-f) 
and Fig. 7(j-l). In Fig. 7(d-f), the Rabi oscillations between |F',0) and IG, 1) cannot be 
observed due to the comparatively large cavity decay rate k. |G, 1) decays and produces a 
photon that escapes the cavity. If we increase the value of the single atom cooperativity, for 
a given number of atoms, the population of the |L, 0) state decreases and the population of 
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FIG. 7. The time variations of the Rabi frequencies of the laser pulses applied to the four-level 
system for N interacting atoms with (a-c) Hi — 2 'kx 72.6 MHz = 1.45^, and (g-i) Hi — 271 x 8 MHz 
= 0.15^. The effective Rabi frequency S{t) is represented by the black solid curve. The red dotted 
curve shows the Rabi frequency of the control laser pulse 12(t) and the blue dashed curve shows the 
cavity coupling strength g. Panels (d-f, j-1) show the corresponding probabilities for the population 
of states |G,0) (black solid), |F^,0) (red dotted), |i?, 0) (blue dashed), |G, 1) (magenta dot-dashed) 
and |T,0) (green dotted) with (d-f) Hi — 1.45^, and (j-1) Hi = 0.15^. Other parameters are the same 
as in Fig. 2. 
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the cavity ground state is enhanced. If we increase the number of atoms, we observe that 
the ground state population is also increased. Thus, the increase in effective cooperativity 
resulting from an increase in the number of atoms plays a constructive role in enhancing the 
efficiency of photon production. In Fig. 7(j-l), the Rabi oscillations between \E, 0) and IG, 1) 
cannot be observed because there is a coherent population transfer from the Rydberg state 
\R, 0) to the cavity ground state |G, 1). The transient population in \E, 0) is approximately 
zero. This reduces the loss due to spontaneous emission into the dummy state |T, 0) and 
improves the photon emission efficiency into the cavity mode [67, 68]. 

We can also investigate the more complicated situation where two Rydberg atoms are 
excited within the blockade radius. The level scheme for this case is shown in Fig. 8. Using 
a sequence of pulses with effective Rabi frequencies Si(t) and S 2 {t) while changing the one 
photon Rydberg level detuning A{t) to compensate the blockade shift Ar, we can excite two 
Rydberg atoms into the state \RR,0) [69, 70]. Then we have two channels: (1) \RR,0) —?• 
\ER, 0) ^ \EE, 0) ^ \E, 1) ^ \G, 2) and (2) \RR, 0) ^ \ER, 0) ^ \R, 1) ^ \E, 1) ^ |G, 2) 
by which |G, 2) can be produced. By applying a control laser pulse with Rabi frequency U(t) 
and the cavity coupling strength g, we can observe emission of two photons from different 
atomic dipoles out of the cavity mode. By using Eq. 16, we can calculate the photon 
correlations. 

In Fig. 9(a, b), we show the Rabi frequencies versus time for the system of = 2 
interacting atoms. Rabi frequency Si{t) drives the first excitation while S 2 {t) drives the 
second excitation. The cavity decay rates are k = lAbg and k = Q.lbg in Fig. 9(a, b) 
respectively. In Fig. 9(c, d), we can see that a 1 Rydberg excitation state is generated by 
using an effective Rabi pulse S\{t). Subsequently, by applying another effective Rabi pulse 
S' 2 (t), we can generate a 2 Rydberg excitation state. This state is similar to a Rydberg 
atom macrodimer [58]. A single off-resonant excitation of the macrodimer can also be used, 
equivalent to a photoassociation of the molecule [71]. By using control laser pulse Q,{t) with 
cavity coupling strength g, we observe that more transient population builds up in \E, 0) and 
\EE,0) compared to |i?, 1) and |G, 1). This situation increases the spontaneous emission 
from |E,0) and \EE,0) to the dummy state |T,0). The photon emission efficiency in the 
weak cavity coupling regime is reduced due to the spontaneous emission. If we increase the 
cooperativity the photon emission efficiency can be improved. Fig. 9(d, f) corresponds to 
the strong coupling regime. In this case, we can populate the 2 Rydberg excitation state and 
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|G,0> |G,1> |G,2> ^ 

(b) Two Rydberg excited atoms 

(a) Level scheme for two Rydberg excitations coupled with cavity 


FIG. 8. (a) Level scheme for two Rydberg excitations with N interacting atoms coupled to a cavity 
with coupling strength g [blue dashed line], S'i(t) [black solid line] and S'2(t) [brown solid line] are 
the effective Rabi frequencies of the laser pulses for one and two Rydberg excitations, respectively. 
fl{t) [red dotted line] is the Rabi frequency of the control laser pulse, (b) Pictorial view of two 
Rydberg excited atoms coupled with the cavity. The excitation is effectively a macrodimer state. 

observe photons emitted from the cavity. We see that the photons are quickly emitted from 
the cavity and photon emission is improved because the transient population of 0) and 
\EE,0) is negligible compare to \R, 1) and [G, 1). The population is coherently transferred 
from the Rydberg state to the cavity ground state, similar to the single excitation case. 

In Fig. 10, for N = 3 interacting atoms, the detunings and blockade shift are the same 
as in Fig. 9. Similar to Fig. 9, the photons are quickly emitted from the cavity and some 
population decays into the dummy state |L,0). The photon emission probability is larger 
in this case because the effective cooperativity is larger. Although small, as we increase the 
number of atoms the photon emission probability is enhanced, due to the effective increase 
in g. 

The expectation value of the dipole-dipole correlation for the case of two Rydberg excited 
atoms is shown in Fig. 11. This quantity is given by 

= {«1<|0+D-|\P> (33) 

and measures the cross-correlation existing between the dipoles belonging to the two atoms 
[64]. Df are the raising and lowering operators for the cavity transition which can be defined 
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FIG. 9. Coherent manipulation of two interacting Rydberg atoms with two excitations, (a, b) The 
time variations of the Rabi frequencies of the laser pulses applied to the four-level system for = 2 
interacting atoms with (a) Hi = 27rx 72.6 MHz = 1.45^, and (b) Hi = 27rx 8 MHz = 0.15^. Rabi 
frequencies Si{t) and 5^2(t) are represented by the black solid and red dotted curves, respectively, 
and generate the first and second Rydberg excitations, respectively. The blue dashed curve shows 
Q{t) and the magenta dot-dashed curve shows g. (c, d) Probability of finding atoms in |G, 0) (black 
solid), |i?, 0) (red dotted), 0) (blue dashed), \ER^0) (magenta dot-dashed) and |iv,0) (green 
dotted), (e, f) Probability of finding atoms in \EE^0) (black line), \E^1) (red line), |G, 2) (blue 
line), |i?, 1) (magenta line), |£’,0) (green line), |G, 1) (brown line). Parameters for the calculation 
are P^ = 27r x 1.4kHz, P^ = 27r x 3 MHz = Po/2, Ac = 0, Aji = 27r x 220MHz, = 27r x 110MHz 
and A(t) changes with time from — 27r x 110 to 55 MHz. 


as D^\g)i = \e)i and D^\e)i = \g)i^ where subscript i represents the atom. In Fig. ll(a,b), 
the normalized dipole-dipole correlation function as a function of time is shown. We see that 
the dipoles are correlated. The laser pulses and cavity coupling time dependencies are the 
same as shown in Fig. 9(a,b). In Fig. 11(a), with spontaneous emission, the prominent dip 
observed is due to the transient population in the cavity excited state that results from 


23 

































































0 1 2 3 4 0 1 2 3 4 1.9 2.2 2.5 

Time (jus) 


FIG. 10. Coherent manipulation of three interacting Rydberg atoms with two excitations, (a, b) 
The time variations of the Rabi frequencies of the laser pulses applied to the four-level system 
for A/^ = 3 interacting atoms with (a) = 27rx 72.6 MHz = 1.45^, and (b) = 27rx 8 MHz 

= 0.15^. Rabi frequencies Si{t) and 5^2(t) are represented by the black solid and red dotted curves, 
respectively, and generate the first and second Rydberg excitations, respectively. The blue dashed 
curve shows and the magenta dot-dashed curve shows g. (c,d) Probability of finding atoms in 
|G, 0) (black solid), |i?, 0) (red dotted), 0) (blue dashed), \ER^0) (magenta dot-dashed) and 
|L,0) (green dotted). (e,f) Probability of finding atoms in \EE^0) (black line), \E^1) (red line), 
|G, 2) (blue line), |i?, 1) (magenta line), \E^0) (green line), |G, 1) (brown line). Other parameters 
are the same as in Fig. 9. 


spontaneous emission outside the cavity. 

Fig. 12(a,b) shows the rate of photon emission efficiency for two interacting Rydberg 
atoms with two excitations for cavity decay rates k. = lA5g and k. = 0.15g respectively. 
In Fig. 12(a,b), the weak cavity coupling case has larger photon emission amplitude than 
the strong coupling regime. With the addition of spontaneous emission, the rate of photon 
emission decreases significantly in the weak cavity coupling case, again due to the transient 
population in the cavity excited state. 
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FIG. 11. Expectation value of the dip ole-dip ole correlation versus time for the case of 2 interacting 
atoms with (a) k = 27tx 72.6 MHz = 1.45^, and (b) k. = 27tx 8 MHz = 0.15^. Other parameters 
are the same as in Fig. 9. 

The second order correlation function of the photons emitted by the two different atoms 
is given by 

where subscripts i,j represent the and atoms respectively. If Log[gj-^^] > 0, radiated 
photons tend to be bnnched while if Log[gj-^^] < 0, the photons are anti-bnnched. This 
function is plotted in Fig. 12(a,b) with and without spontaneous emission as a function of 
time for k = lA5g and k = O.lbg. In these cases, we observe that Log[gj-^^] is smaller than 
zero for short times when spontaneous emission is neglected. This indicates that the photons 
are anti-bunched at short times when the photons are preferentially emitted. The value of 
Log[g|^^] is greater than zero for longer times, indicating the emission of photons is bunched 
after the initial photon pnlse. If we consider the spontaneous emission into the dnmmy 
state |T, 0) then Log[g^^^] is always greater than zero. Log[gj-^^] > 0 implies the emission of 
photons is always bunched for the more realistic case. The additional dephasing destroys 
the interferences of the 2 oscillating dipoles in this case. 

To dehne the entanglement and purity of the quantum systems, there are several mea¬ 
sures. One of them is the von Neumann entropy of the quantum system. This is a very 
sensitive measnre of the purity of the quantum system and it determines the degree of cor¬ 
relation between the collective excitations and photon states [72-74]. The von Neumann 
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FIG. 12. (a,b) The rate of photon emission efficiency for 2 interacting atoms with two Rydberg 
excitations. (c,d) The logarithm of the second order correlation function Log[g^ 2 h of photons 
emitted by two different atoms with two Rydberg excitations versus time. (e,f) Variation of total 
entropy of the system with time. (a,c,e) k . = 27rx 72.6 MHz = 1.45^ and (b,d,f) k . = 27rx 8 MHz 
= 0.15^. Other parameters are the same as in Fig. 9. 


entropy of the atoms-cavity system can be defined as 


with density operator 


j 


P= - 

j 


(35) 


(36) 
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where the r/j^s are the diagonalized eigenvalues of the system and j is the number of states in 
the Hilbert space. Since our system is multi-partite, and there is no general measure to study 
the entanglement for multi-partite system, we approximate our system as a bipartite system 
of atoms and photons. Then, the sum of the entropies of the 2 sub-systems atoms, S{pa), and 
photons, S{pp), should be larger than the total entropy of the system, S{pa) + S{pp) > S{p). 
The bound for the entropy for these two interacting sub-systems can be dehned by the Araki 
and Lieb [75] triangle inequality, 

|S(ft) - ■S'(Pp)l < S{p) < |S(p„) + S(p,)|. (37) 

The total entropy of the system is plotted in Fig. 12 taking into account and also ne¬ 
glecting spontaneous emission as a function of time. Fig. 12(e,f) shows the variation of total 
entropy with cavity decay rates k = lAhg and k, = 0.15^'. As photons are emitted from 
the cavity, the entropy of the system first increases but then decreases due to the loss of 
the photons. In Fig. 12(e,f), the entropy of the system is larger for the realistic case where 
spontaneous emission out of the cavity is considered. In both cases, after the complete 
emission of the photons, the entropy of the photonic sub-system has reached zero with time 
and system is in a pure entangled atomic state [76, 77]. 


VII. DISCUSSION 

In the literature, there are several different ways to generate single photons [42, 43, 78- 
80]. One way is to use four-wave mixing [50]. The approach we have investigated in this 
work is superficially similar to the four-wave mixing method because of the atomic level 
scheme. One fundamental difference is that the method described here does not require 
phase matching. In the four-wave mixing method, the phase matching condition needs to 
be fulfilled to observe the photon which will be emitted in a diffraction-limited solid angle. 
The approach described in this paper uses an alternative way to generate single photons by 
preparing a single collective state using Rydberg atom blockade in a cavity. The lasers are 
time dependent in such a way that the preparation, and control are not present at the same 
time. Technologically, Rydberg blockade allows one to deterministically control the number 
of excitations in the cavity and the cross-section for coupling to the cavity is enhanced 
by the collective nature of the initial state. This addresses some of the problems inherent 
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in trapping single atoms in an optical cavity to generate single photons. The diamond 
configuration enables one to isolate the quantum state preparation from the single photon 
output. The output is in a single mode determined by the properties of the optical cavity. 

A similar atomic level scheme has been studied for collective excitation for Raman and 
superradiant lasers [39-41]. In these systems, N atoms are excited collectively without con¬ 
sidering atom-atom interactions. Our model shares the structure of the excitation scheme, 
but Rydberg blockade prepares an N atom entangled state with a controlled number of ex¬ 
citations or controlled number of singly excited superatoms. In a more general case, spatial 
correlations can be written into the sample. This is quite different from these laser schemes 
as our approach produces a state similar to a micromaser but with deterministic control 
over the number of excitations in the cavity. The system we described can be useful for 
studying superradiant behavior because a collective state can be deterministically prepared 
as the initial state. The cavity coupling can be used to analyze the collective quantum state 
by detecting the emitted radiation [64]. 

Most of the models discussed previously for Rydberg atoms in a cavity use ground state 
atoms as an initial state. The JC model for the atoms initially in the ground state in an 
optical cavity with a single Rydberg excitation as an essential piece of the dynamics has 
been proposed in Ref. [38]. This model is different from our approach. In our model, we 
first excite a single collective state by using the Rydberg blockade mechanism and then 
by manipulating a control laser field, generate photons from the cavity. In addition, the 
preparation lasers are all far detuned from the cavity field frequency. Now, if we consider 
the Rydberg excited state as an initial state of the collective atoms-cavity system, then the 
photon emission efficiency can be increased due to the collectively enhanced cooperativity 
and can be further improved once the strong coupling regime is achieved due to a negligible 
transient population in the cavity excited state [64]. 

We have neglected several dephasing processes in our analysis because we anticipate these 
effects are small, or can be made small with a judicious choice of parameters. The effects that 
are of most concern are the spatial dependence of the atoms-cavity coupling, dependence 
of the Rydberg atom interactions on internuclear separation, motion of the atoms, and 
interaction between Rydberg excitations when there are 2 excitations simultaneously in the 
cavity. Taking these effects into account with the straightforward Liouville equations we 
have used in our calculations is involved, so we estimate their effect on our current results. 
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Light pressure forces on the atoms in the cavity are negligible because we have assumed 
the cavity is on resonance with the 5^5'i/2 ^ 5^P3/2 transition for times after the atoms 
are excited to the h‘^Pz /2 state, the exception being Fig. 3 where we show a detuned cavity 
to demonstrate that our solution is producing expected results. There is no asymmetry in 
the absorption of the photons in the cavity under these conditions and the photon recoil is 
small. The recoil temperature of ®^Rb is 362 nK. Additionally, for a large fraction of the 
time the atoms are interacting with the cavity, there are no photons inside it. 

We have run our calculations with additional dephasing rates by coupling the singly 
and doubly excited Rydberg states to the dummy state. The results of these calculations 
with additional dephasing provide an estimate of the effects of dephasing on the collective 
excitations. In these calculations, we used the same rates for each transition. Both processes 
were modeled with the Lindblad operator. We examined additional dephasing rates up to 
50 kHz. The results for the photon production probability are reduced by around ~ 4% for 
50 kHz of dephasing on each transition for the excitation timescales used in this work in the 
2 excitation case. The reduction in the overall photon production scales linearly as these 
dephasing rates are increased. 

Practically, ultracold samples at temperatures ~ 100 nK can be generated and used for 
an experiment to reduce the decoherence resulting from thermal motion, although photon 
recoil can limit the temperature to 362 nK. At temperatures of ~ 1 /uK, the atoms only 
move ~ 100 nm in 10 /rs. Because this is a small fraction of a blockade radius, free evolution 
of the atoms can be neglected at ~ 1 /uK. The corresponding Doppler shifts of the atoms 
at ~ 1/iK are ~ 27 r x 10 kHz. Doppler shifts on this order can lead to reductions in the 
photon production rate of ~ 1%. The laser spectral band widths and the inverse of the times 
over which the collective states freely evolve are broader than the Doppler widths of the 
atoms that make up the collective excitations. In our case, for very low, but experimentally 
achievable temperatures, the frozen gas approximation is valid to around 1% in the calculated 
photon production probability. 

In general, multiple excitation states within the blockade volume and states consisting 
of more than 1 superatom will have spatial correlations [81-85] because the interaction 
between the Rydberg atoms depends on their internuclear separation and the orientation 
of the internuclear axis relative to any helds. However, as mentioned earlier in the paper, 
2 Rydberg atoms can be excited around an avoided crossing where there exists a broad 
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region where the energy has weak dependence on the internuclear separation. These types 
of potentials can be observed when the multi-level nature of the Rydberg atom interactions 
and external magnetic and electric fields are considered. In fact, the fiat regions of the 
Rydberg potentials, so called ‘stationary regions’, are regions where Rydberg atom pair 
excitations occur most readily. The acceleration of the atoms after excitation is determined 
by the gradient of the interaction potential and has been observed experimentally [2, 58, 59]. 
The acceleration on these types of interaction potentials can be negligible as demonstrated in 
these works - in these experiments, we had to wait ~ 200 /is to observe a difference between a 
thermally expanding cloud of cold atoms at 80 /rK and the Rydberg atom collision partners 
[59]. Scaling this time to a 1/rK temperature gives a time of ~ 20/rs. Similarly, if we 
assume the atoms constantly accelerate to the final velocity measured in [58, 59], 17 cms“^, 
in ~ 100 /rs, it takes ~ 5 ;us to reach a thermal velocity that corresponds to 1 /rK. The time 
over which the 2 Rydberg excitations exist in our case is only ~ 1 fis, so this effect can 
be negligible with an optimized choice of interacting Rydberg state, or use of electric and 
magnetic fields to engineer the interaction potentials [58, 86, 87]. Recall that the acceleration 
will cease when the 2 Rydberg excitation state is projected onto the non-interacting atom 
collective state. For interaction potentials where there is significant radial dependence, the 
efficiency of the 2 excitation dynamics will be reduced and there will be spatial correlations 
in the state that is excited. Changes in the efficiency can depend strongly on the Rydberg 
atom interaction potentials. 

For the calculations in this paper on 2 excitation states, the collective excitations are 
not independent of one another but are analogous to a doubly excited Rydberg molecule 
or macrodimer. In this case, the dephasing rates of these entities should be approximately 
equal to those occurring in a 2 independent excitation case where 2 superatoms are present, 
provided there are no additional open molecular decay channels. In many cases, the decay 
time for the macrodimer is determined by the radiative decay of the Rydberg atoms that 
make up the molecule [71, 87]. Ground state-Rydberg atom molecules are not important 
because the internuclear spacings are greater than the size of the Rydberg orbital [88, 89]. 

The cavity coupling strength also has spatial dependence which must be considered. 
The size of the collective excitation is very large compared to the wavelength of the cavity 
transition. The relative sizes can be quantified by comparing the blockade radius for n = 90 
to the wavelength corresponding to the cavity resonance, 9/im ^ 780 nm. The cavity 
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coupling strength is effectively averaged over the spatial dependence of the cavity field 
mode, reducing the average cavity coupling strength by \/2 in this simple argument. A 
rough estimate places half the atoms at nodes and the other half at antinodes. Using this 
approximation to estimate the spread of atomic energy shifts for the atom-cavity coupling, 
yields a spread of atom energies of ~ \fNg. It is interesting to explore the dephasing caused 
by the interaction of a collective excitation that is spatially larger than the modulation 
in the cavity coupling in more detail, but such an investigation is outside the scope of the 
work here. Numerically, this can be approached using quantum Monte Carlo methods which 
are computationally intensive but possible [38]. For a collective state spread over a large 
number cavity resonant wavelengths, the overall effect is primarily to reduce the number 
of atoms participating in the atoms-cavity dynamics since atoms in the antinodes of the 
cavity are affected by the cavity coupling while those at nodes essentially behave as if in 
free-space. The atoms located at antinodes dominate the dynamics because they are coupled 
more strongly to the cavity. The spatial dependence will lead to fluctuations of the atoms- 
cavity coupling parameters that will decrease as the number of atoms in the collective state 
increases. Although our calculations are for low atom number collective states, the results 
scale to higher atom number collective states, at least with respect to the spatial dependence 
of the cavity mode. 

The discussion of the dephasing mechanisms here points to a common factor that is 
important for a broad range of Rydberg atom quantum optics experiments. The product of 
the time the atoms exist in a particular collective state and the dephasing rate due to atomic 
motion, including heating mechanisms that arise in an experimental scheme, needs to be 
small. As a consequence, our method requires cold samples, ~ 1 /rK, and fast excitation 
dynamics, < 1 /is. In addition, to create multiple excitation states, we require potentials 
that are engineered to reduce the acceleration of atom pairs after excitation and provide 
a broad range of internuclear separations over which atom pairs can be excited, if spatial 
correlations are undesirable. Speeding up the excitation processes used in our scheme can 
compensate for some of the losses due to the additional dephasing mechanisms if, as in the 
case investigated, the Rydberg atom interaction strength is large enough to accommodate 
an increase in the Rabi frequencies of the laser pulses. 

The Purcell factors and cooperativities we have used can be justified on the grounds that 
some experiments have achieved cavities with sufficient finesse and mode volumes to at least 
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approach, or realize, the parameter values {Fp, g, k) used here [90-92], Lower cooperativity 
factors, of course, will reduce the elRciency of the approach that we have presented, however, 
a smaller cooperativity factor does not preclude the qualitative observation of the proposed 
scheme to deterministically synthesizing quantum states in a resonant, high-finesse cavity. 
Realizing this goal is interesting and can motivate improvements in the cavity design. 


To summarize, we showed that Rydberg blockade allows one to deterministically prepare 
a fixed number of collective excitations inside a cavity using our model. This is impos¬ 
sible for the case of direct excitation of atoms inside a cavity, since, for noninteracting 
atoms, the excitation is probabilistic. In the case of 2 excitation states generated by the 
blockade mechanism, we can deterministically observe interesting dynamics and analyze the 
correlations between two, or more, atoms and the photons emitted from the cavity. The 2 
excitation states we examined in this paper can be understood as macrodimers excited on 
flat relatively featureless interatomic potentials that are subsequently coherently transferred 
to noninteracting, collective excitation states. Both the excitation and transfer processes 
produce collective states as the excitations are shared by the atoms in the sample. Using 
macrodimer states that are not featureless can be used to introduce spatial correlations into 
the collective excitations. A detailed treatment of the spatial correlations introduced by 
different interatomic potentials lies outside the scope of this paper, but is a future goal. The 
engineering of spatial correlations by using Rydberg atom interactions in a collective en¬ 
semble is interesting and could be useful for experiments in quantum optics, for example in 
studying superradiance. In this paper, we showed that the collective state dynamics can be 
controlled deterministically by demonstrating results for collective dynamics that are similar 
to the TC model for two excited atoms. We have observed anti-bunched photon emission 
when we neglect spontaneous emission out of the cavity. With spontaneous emission the 
emitted fight is bunched in the weak {k > g > T±) and strong (g > k > r±) coupling regime. 
Similarly, the entropy of the system decreases with time, while the second order correlation 
function of the photons emitted by the 2 collective excitations is enhanced. The efficiency 
of photon production is enhanced by increasing the effective cooperativity by increasing g 
through the cavity parameters or the number of atoms, as would be expected. 
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VIII. CONCLUSION 


We have investigated a four-level system of N entangled atoms placed inside a high-finesse 
cavity. The Rydberg blockade mechanism plays an important role in preparing the atoms 
in a collective state. We have used numerical simulations to show that one can deterministi¬ 
cally excite collective N atom superposition states using Rydberg blockade and laser pulses. 
Rydberg blockade allows one to deterministically prepare a fixed number of collective exci¬ 
tations inside a cavity. This is impossible for the case of direct excitation of atoms inside a 
cavity since that process is probabilistic. In the case of 2 collective excitation states gener¬ 
ated by the blockade mechanism, we can deterministically observe bunched photon emission 
from the cavity similar to the case of 2 excited atoms initially prepared in a cavity. The 
scheme is useful for the generation of single photons as well as other types of quantum light 
fields. It could also be useful for coherent optical manipulation in quantum information pro¬ 
cessing and the study of interesting quantum optical phenomena such as superradiance [93]. 
The rates used for the calculations are realistic and experimentally achievable. The most 
important dephasing mechanisms were included or addressed. This system merits further 
theoretical and experimental study because it demonstrates flexibility for deterministically 
synthesizing collective quantum states that can be analyzed by interrogating the photons 
leaving the cavity. 
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